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Abs tr act 

Atmospheric extinction in wideband photometry Is examined 
both analytically and through numerical simulations. If 
the derivatives that appear In the Strttmgren-K Ing theory 
are estimated carefully# it appears that wideband 
measurements can be transformed to outside the atmosphere 
with errors no greater than a m i 1 1 i magn I tude • A numerical- 
analysis approach is used to estimate derivatives of both 
the stellar and atmospheric-extinction spectra# avoiding 
previous assumptions that tne extinction follows a power 
law. However# it is essential to satisfy the requirements 
of the sampling theorem# to keep aliasing errors small. 
Typically# this means tnat band separations cannot exceed 
half of the full width at half-peak response. Further work 
is needed to examine higher-order effects# which may well 
be significant. 
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I ntr oduc ti on 


The reduction of wideband photometric observations to 
outside the atmosphere# and thence to some standard system# 
are Incompletely solved proDlems of long standing* As Mas 
shoMn by Bengt Strdmgren [19373# and emphasized by Ivan King 
C19523# these problems are Intimately related; for He may 
regard the extinction correction as a color transformation 
that depends on sir mass* 

This problem is hardly new* In the very earliest 
photoelectric photometry# GuthnicK and Prager [19143 found 
that "The correction for extinction is one of the most 
difficult problems for such exact measurements as can be 
reached by photoelectric methods.... The extinction is 
dependent on the spectral type to a high degree. Under 
normal transparency conditions# the ratio of the 
photoelectric to the visual extinction is about 2.2 for the 
middle of class B# about 2.0 for class A# etc.# about 1.3 
for class Ma. It Mill apparently turn out later that these 
factors themselves are also functions of the zenith 
d i st a nee • • • •" 

Indeed# Forbes [18423 had already found "That the 
tendency to absorption through Increasing thicknesses 
of air is a diminishing one.... Hence the amount of vertical 
transmission nas alMays hitherto been greatly overrated# 


216 



or the value of extra atmospheric solar radiation greatly 
unde r r at ed • • • • The physical cause of this I a* of absorption 
appears to oe the non-homogeneity of the incident rays****” 

The first analytical treatment of these problems Mas 
made by ri.S tr Bmgren [1937J; Mho showed that if the spectral 
sensitivity curve of a photometric instrument is fairly 
narrow; a Taylor series expansion of the stellar spectral 
irradiance curve about the instrumental centroid 
wavelength allows the extinction to be expressed In terms 
of the mon ocnr om at i c extinction at this wavelength; with 
a correction term proportional to tne square of the 
instrumental bandwidth. This approach was developed 
further by King [1952]; using a more compact notation. 

King's paper Is required reading for anyone who wants to 
understand he ter ochromat i c extinction. 

The correction terms involve first and second derivatives 
of both the stellar spectral irradiance and (in the 
extinction problem) tne atmospneric transmission. As is 
snown by King ti.952] and Young [1974]; the first derivatives 
may oe approximated by the colors of the stars and the 
reddening of tne atmosphere; respectively. King [1952] 
shows that tne second derivative for the stars can be 
neglected; and that — under certain assumptions tnat Mill 
be discussed below -- tne second derivative for the 
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atmosphere can be related to the first derivative. 


Although the corrections are traditionally made by using 

color indices as independent variables* Cousins and Jones 

Q9763 found that "no equation involving 8-V and U-B only 

Mill predict the extinction correctly for all luminosity 

types and different degrees of r eddeni ng. • * • The difference 
m 

can exceed o.Oi. Without more information* direct or 
Inferred* no rigorous colour correction is possible either 
for extinction or for colour tr ans f ormati on. • • •" Similar 
conclusions Mere reached independently by Mandwewala [19763; 
see also Blanco [19i>73* and hig.14 of Young [19743. 

Some years ago* the accepted Mlsdom was that these 
difficulties were due to the great width of the UBV passbands* 
and that I ntermedi ate-passba nd systems such as uvby would 
prove far superior. However* bandwidth effects are 
proportional to the square of the passband width* which 
is about 3 times narrower for uvby than for UBV* hence* if 
this were the only problem* such difficulties should be 
nearly an order of magnitude smaller for uvby than for UBV. 

But* after the most strenuous efforts at 
standardization* Olsen [19833 found typical systematic 
differences oetween 4-color data from northern and southern 
stations on the order of 0.004 mag.* and unexplained systematic 
errors of several nundredtns of a magnitude for a number 
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of individual stars. These figures are only about a factor 
of 2 better than has been done in comparably careful UBV 
photometry, lianfroid [1985] says that "reduction of many 
observing runs in the uvby system with various equipment 
shows that errors as nigh as .05 magnitude# and more# are 
not uncommon." Indeed# Manfroid and Sterken [19871 have 
recently shown systematic errors as large as a third of a 
magnitude in careful uvby observations# taken at a good 
site (La Sillai# calibrated with dozens of standard stars# 
and reduced by reliable techniques. 

Furthermore# the most precise published photometry 
appears to be tnat done In the Geneva system Csee Fig.2 of 
Young# 1964a]# which has been reproduced to better than 
0.003 mag. for we I l-ooser ved stars Ccf. Table IV of Rufener# 
19811. This precision is all the more remarkable because 
the Geneva bands are comparable in width to those of the UBV 
system# and# liKe those of UBV# are defined by glass rather 
than interference filters and Dy the tail of the 
photocathode response function — unlike the supposedly 
superior and " f i I ter ~def i ne d" uvoy bands. Clearly# there 
Is more to precision than bandwidth alone. 

Because the Geneva worKers have been extremely careful 
to determine ana use derivatives correctly# one suspects 
that a careful examination of the derivative problems 
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would be helpful. 


I have already suggested [Young# 197*; 


198*b] that one problem with the traditional approach is the 
relatively poor approximation to the first derivatives 
provided by spectrally unoersampled data; unfortunately# 
all existing pnotometrlc systems violate the r equ i r ements 
Imposed by the sampling theorem (though the Geneva system 
is less unoersampled than most). 

This has been confirmed by Manfroid [19851# who shows 
that precision can be greatly improved by using a second 
set of filters slightly displaced from the normal set# to 
estimate more accurate first derivatives. I show below 
that both first and second derivatives can be determined 
accurately# and that such improved estimates lead to much 
more accurate extinction and transformation corrections. 

In the past# precision on the order of 1 percent was 
regarded as "good enough" for most problems# though this 
involved an element of circular reasonings the I acK of 
better measurements prevented anyone from even considering 
investigations that required much better than 1 percent 
precision. Today# however# there are a number of problems 
whose photometric study clearly requires precision on 
the order of one m i I I i magni t ode or better: the detection 
of planetary systems tborucKi# 19a*]; stellar seismology 
[Fossat# 198*; Hudson# 198*]; inventorying the Sun*s comet 


220 



cloud [Meinel and Meinel# 1986]; nonlinear dynamics of 
pulsating white dwarfs [Auvergne and Baglin# 19863* We may 
take tnis as a nominal goal to reach — roughly an order of 
magnitude oetter than current practice* Heintze et al* 

C1984] and Schml dt-Kal er [1984] have even suggested that 
still smaller errors could be obtained from the ground* 

It is generally accepted that major advances in ground- 
based photometry will require multi-channel techniques# 
to remove atmospheric transparency variations* However# 
the price that must oe paid is the problem of calibrating 
the different channels against one another# this is 
essentially the transformation problem* whether we use 
multichannel instrumentation or not# we cannot expect to 
do enormously better than 1% if the model used to represent 
tne heterochromatic extinction is no better than 11. 

A number of distinct issues must be resolved* The question 
of adequate spectral sampling was raised above* But# even 
with properly sampled data# how should the derivatives be 
extracted from the data? And# even before these questions 
can be discussed# there is a conceptual problem with the 
derivatives that needs clarification# so let us begin with 
It. 
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Meaning of the Derivatives 


Outline of tne fneory 

To see wny there is a problem with these derivatives# let 
us review where they come from. Let I(X) be the stellar 
spectral irradiance function; let t(X# z) be the atmospheric 
transmission function 

t ( X # z) - exp[-A(X ) H(z)/1.08bJ # (1) 

which King rather inconveniently calls Q(X). Here 

A( X ) is the wavelength-dependent extinction coefficient 
In magnitudes per air mass; rt(z> Is the air-mass function of 
zenith distance z# discussed at length by Young C1974]# and 
1,08b Is short for 2.b/ln(10) * 1.085 7362...# which is the 
conversion factor oetween natural logarithms and logs to 
the base Z.olL (i.e.# magnitudes). As the following 
discussion focuses on the wavelength dependence# we shall 
usually omit the z dependence of t(X). 

If the response function of the Instrument (including 
the telescope optics) is k( X )# the quantity measured when 
we observe tnls star is 

L ■ J U X ) t( X ) k( \ ) d X • <2 ) 

King [1952J splits the integrand into an instrumental part# 

K( X ) # and the rest# 

i( \ ) * I ( X ) t( X )# (3) 

nhich changes from one observation to the next* 
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The heart of the S trOm^r en-K log method is to expand S(A) 

in a Taylor series about some central wavelength# X * 

o 

s(X> * s ( X ) *smX > (X - X > ♦ - s"c X ) (X-X ) 2 ♦ ...» 

o o o 2 o o 


where primes denote wavelength derivatives. The integral 
(2) can then be done termwisei 

L * S J R ( X ) dX + $' j < X ~ X ^ ) R < )s > dX 
+ - ( X - X > 2 *< X ) d X + ... . 

I J o 


(5) 


Thus# the part of the measured light that changes is 

expressed in terms of the transmitted spectral distribution 

S and its derivatives# evaluated at X # and the invariant 

o 

instrumental influence is expressed by the moments of the 

response function k about X • how# If X is chosen to be 

o o 


the centroid wavelength 

f X R( X ) d X 

T *<X~ , " d x 


X 


(6) 


the S' term in (b) vanishes# and we have only 

L * / Kl X > dX [S + 2 S"J * 


(7) 


where JUL is tne normalized second central moment of fc(X )* 

' 2 


>U- 


f ( X - X ) 2 k{ X ) d X 

c J o 

! Tr< X ) d X 


( 6 ) 
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King tnen expresses S H in terms of the derivatives 


of I ( X ) and t( X ) » 

* 1" t ♦ 2 I • t « ♦ it" # 


( 9 ) 


which allows the measurement to be expressed in magnitudes 
as f o 1 1 o ws i 

m * m ♦ A( X )rt(z) - i . 086 In (i ♦ x) # (10) 

ob s o o 

where m is the observed instrumental magnitude# m is 
obs o 

the magnitude that would have been measured outside the 
atmosphere# and 


I <*) I 


V 


I H /I ) - 


1.086 


2 {‘V 


a"/a) ♦ 


X 


( X A'/a) l 2( X l'/n - — -( X A ' / A ) ] \T| . (11) 

1.08O /il 


All the parenthetical expressions are evaluated at X • 

o 

The Derivative Problem 


Now# the difficulty is to assign proper significance to 

the derivatives 1* and I M # evaluated at X • What is intended 

o 

is obviously not what is sai a literally# for stellar spectra 

are cluttered with absorption lines. If X happens to 

o 

fall on the sloping side of a line# I 1 will be enormous# and 
obviously untypical of the general run of the spectrum in 
the region a few honored Angstroms wide that is of interest 
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In photometry, doth Strflmgren and King neatly sidestepped 
this question by considering only blacK bodies in their 


examples. Yet we must deal with real stars# whose spectra 
plainly cannot be represented accurately by any low-order 
Taylor series. 

Evidently# the only sensible physical interpretation 
of these derivatives is that they refer not to the true 
stellar energy distribution# which fluctuates wildly# but 
to some smooth function that approximates It# and that has 
well-behaved derivatives. As only first and second 
derivatives appear in 1X1)# we may suppose that a parabolic 
approximation is used. 

if we write the true stellar spectral irradiance as 

it X ) « i ( X ) + I ( X ) # (12) 

s r 

where I and I are the smoothed intensity and the 
s r 

remainder# after subtracting the smoothed irradiance from 
the true one# we want to choose the smoothing so that all 
the preceding equations are true when 1^ and its 
derivatives are used in place of the true I ( X ) # I * # and I". 
In particular# suppose we use (121 in (2)# so that 

L * J I ( X ) t C X ) k ( X ) d \ + Jl f (X) 1 1 X ) R ( X > d X . 

Then the second term in (13) must be zero# if King’s 

formulae are to be true for I • To make 

s 


( 13) 
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■ 0 # 


(14) 


J 1 ( X ) t( X ) R(X ) d X 

~ / r 

expand t( X ) in a Taylor series about X [chosen according 


o 

to (6)1. This expansion Is generally well-defined# because 
t< X ) is quite smooth# provided that we avoid spectral regions 
containing sharp molecular absorption lines. Then 


t()l ) j I i X ) R C X I d X ♦ t*(X ) j(X - X ) I (X) R ( X ) dX ♦ 

o ^ r o or 

\ f(A * < X " X > 2 i(X)R(X)d\ - o. 

J A I A • * ' 


0 J 


o r 


(15) 


Because t( ) and its derivatives depend on air mass# and 
o 

change from night to night# (15) can be generally true only 

if each of these integrals vanishes. For a given star# and 

hence a given I(X )# equating each integral to zero provides 

three conditions on 1 # and hence on 1 . These are the 

r s 

three conditions required to choose the parabolic function 
1^ ( \ ) uni que I y . 


In what follows# 1 assume that this choice has been 

made# so that all the formulae refer to 1 and its well- 

s 

defined derivatives# instead of to 1. The subscript s will 
be suppressed# but must be understood to be present 
throu ghout • 

Evaluating the Derivatives 

Although the discussion above clarifies the meaning of 
the derivatives that appear In (11)# it does not provide any 
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way to evaluate them in practice* The whole point of 
broadband photometry is to avoid spreading the light of 
faint stars out into spectra# which would be necessary to 
evaluate the moment integrals in (15). We need to obtain 
the derivatives in (11) from the photometric data themselves. 

King [19621 used a color index to estimate I'# and this has 
been the traditional practice ever since. He pointed out 
that I" acts primarily as a zero-point shift that Is 
Independent of air mass# so that for many purposes It can 
be ignored. Young and Irvine [19671 and Young [197A1 use 
the reddening power of the atmosphere to estimate A*> this 
Is exactly analogous to using a stellar color index to 
estl m ate 1 1 • 

However# the A" term cannot be ignored. King [19621 argued 

that the w a ve I eng th— dependen t part of the extinction is 

\ -4 

mainly Rayleigh scattering# and hence proportional to h . 
Because A' appears in (11) only In the combination ( X A’/A )# 
which is tne logarithmic derivative of A ( X ) » and hence 
a constant if a(A) is a power law# King argued that the value 
of this constant is -4. His assumption of power-law extinction 
also allowed the second derivative A" to be expressed in 
terms of A ' . 

If we consider the scattering part of the extinction# 
even the Rayleigh scattering is not exactly a power law# 
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because the dispersions of the refractivity and the 
anisotropy of the po i ar i zab t i i ty inane the molecular- 
scattering extinction steeper than A * (Young# 1982} 

Bates# 19b4# Nicoiet# 1984]* On the other nand# the aerosol 
extinction is much flatter# and is usually near 1/A • 

Because the aerosol extinction dominates at long 
navel engths# and the molecular scattering at short 
wavelengths# the logarithmic derivative of A ( X ) is closer 
to -4 In the violet and closer to -1 in the red. 

On top of this# there is very strong absorption by ozone 
below about 350 nm# and more than If. absorption In the 
Chappuis bands between about 500 and 680 nm (Vigroux# 1953} 
Inn and Tanaka# 1953# Griggs# 1966 J. This band absorbs more 
than 0.1 percent in the zenith between about 450 and 650 nm. 

Thus# A(A ) cannot be regarded as a power law function 
In accurate worn# despite tne pedagogical utility of this 
crude approximation. However# If A( X ) is not a power 
law# then not only is (X A'/A) wavelength-dependent# but 
we cannot express the second derivative in terms of the 
first. In fact# the A M term 


X 


2 

A 


A 


2 

d (In A ) 
d( In A ) 2 



(16) 


involves both the first and the second logarithmic 
derivatives of A ( A- )• If A( A ) were a power law# Its second 
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logarithmic derivative would vanish* and (16) would provide 
a very simple estimate for the A N term in (11)* In reality* 
we must estimate the second logarithmic derivative of A( A ) 
in addition to the first* 

In 19 5<d* only one color index (the old international 
color index) was in common use* and the plethora of modern 
multicolor systems had not yet been invented* Even the UBV 
system had not been formally introduced* Thus it was quite 
natural for King to try to make a single color index do 
everything* only one was available* But today* we need not 
be so restricted* 

If we measure only two points on a function* we can 
fit a straight line through them* This linear fit allows 
us to estimate ootn the function and its slope at any point* 
If we have ‘three data* we can fit a parabola* and determine 
the curvature (i*e** a second derivative) as well* Even if 
the data are unequally spaced* standard techniques of 
numerical analysis allow us to find these derivatives at 
any point* 4 

I propose to adopt this nume r i ca I — ana I y s i s point of view* 
so as to estimate the A" term* Though this may appear novel* 
it is really quite similar to what is already done for 
stars in a number of multicolor systems* If the bands are 
nearly equally spacea in wavelength* we can use one color 
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Index to estimate the slope of a star's spectrum# and the 


difference of two neighboring color indices to estimate 


its curvature* In fact# such curvature indices as m^ and 
c in the uvby system are already quite familiar to 


1 

photometr i sts. 

If three bands are exactly equally spaced in wavelength# 
or three samples of any function are equally spaced in its 
argument# it is well known that the first derivative at 
the central point is better estimated by the slope between 
the two end points than by the slope between the central 
one and eitner of its neighoors# provided that the samples 
are sufficiently close together tnat higher-order terms 
can be neglected* (This fact is used in the reduction of 
photometry in the Geneva system# this must surely be another 
reason for tne excellent precision of tire published Geneva 
results*) The reason is simply that three points allow 
parabolic approximation# which is generally better than 
the linear approximation through two neighboring points* 

If the tnree points are unequally spaced# the derivative 
of the function can still be expressed as a simple weighted 
sum of the tnree ordinates# following standard Lagrangian 
interpolation methods. In the general case# the weight of 
the central point is not zero# as it is for equal spacing. 
The details of the derivation are given in tne Appendix* 
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However# this process only makes sense if the bands 
overlap enough that a polynomial passing through their 
average intensities at their effective wavelengths is a 
good representation of the smoothed spectral irradlance 
function# as defined above. In other words# the smoothed 
irradiance must be sampled at intervals (i.e.# band spacings) 
close enough to satisfy the sampling theorem [Young# 19741. 

No existing pnotometric system does this# though the Geneva 
system comes close. 

Numerical Simulations 

To illustrate the improvements possible with adequate 
sampling and accurate data reduction# 1 have done several 
simple numerical simulations of wideband photometry and 
reductions. For simplicity# I used symmetrical passbands# 
to Keep third-order effects negligible. Both inherently 
smooth spectra (a set of black bodies# and a set of artificial 
cont I nua parabolic in the logarithm of spectral irradlance) 
and realistic spectra taken from the tables of Gunn and 
Stryker C1983J were multiplied Dy standard atmospheric 
transmission functions for 1.0# 1.5# 2.0# and 2.5 air 
masses; multiplied by cosine-squared response functions 
500 Angstroms wide at naif maximum; and integrated# to give 
synthetic observational data. Similar calculations without 
the atmospheric transmission gave true extra-atmospheric 
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values for eacn Instrumental system 


In eacn case# tne central band of a 3-band system was 
kept fixed at 450o Angstroms; band spacings of 100# 200# 

300# and 500 Angstroms xere used for the stars. The columns of 
Table I show the standard deviation (root— mean— square 
residual per degree of freedom); the maximum residual in 
the middle (4500 A) band — — a rough measure of internal error; 
and the maximum error in the extra-atmospheric magnitude 
in the middle oand calculated from the fitted parameters. 

This last column is a rough estimate of external error. 

The taole concentrates on the results for the middle 
band# because it is the same for all cases; the outer bands 
move as the band spacing changes# and so are not strictly 
comparable from case to case. Nevertheless# It is worth 
remarking that the errors in the s hor test-wa ve I eng t h band 
are about double tnose for the middle band so long as the 
spectra are smooth# or the bands are closer than 300 Angstroms. 
For 300 A spacing# the errors at the shortest wavelength are 
about 4 times tnose for the middle band; for 500 A spacing# 
they are aoout 5 times larger. 

Thus# the Table suggests that, sampling for real stellar 
spectra is adequate at 200 A spacing; marginal at 300 A; and 
wholly inadequate at 50o A# wnicn Is the full width at half- 
maximum of the bands. The aliasing errors are thus small 
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for band spacings below half the FkHM/ and Increase rapidly 
at larger spacings* Furthermore/ as only a dozen stars have 
been used# and we may expect the worst error from a larger 
set of stars to be somewhat larger/ It appears that 
ml 1 1 I magni tude accuracy can be achieved with bands as 
broad as 500 A if their spacing is about 250 A# 

The reader should bear in mind that these bands were 
perfectly symmetrical/ so that third— order terms (which 
involve the third central moment of the passband) were 
eliminated. Real filters always produce markedly 
asymmetrical passbands/ so we may well need to include the 
next-order terms in the expansions* Unfortunately/ the 
resources available for this work did not allow a thorough 
investigation of the spacing required to reach a given 
level of precision with realizable passbands. 

As the atmosphere is part of the Instrumental system/ 
the success in transforming tnese pseudo-oos e r vat I ons 
from inside to outside the atmosphere to m 1 I I 1 mag n I tude 
accuracy suggests that the transformation problem between 
different instrumental systems can also be satisfied with 
properly sampleo data/ using this same numerical- 
interpolation approach. 

C one I usl ons 

Bandwidth effects do not seem to be a serious I imitation 
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to the precision and accuracy of broadband photometry# if 
they are modelled correctly. This requires both a more 
detailed understanding of the derivatives that appear in 
the classical :>tr9mgren-K ing type of analysis# and the use 
of ue I l-estao I i shed numerical-analysis methods to determine 
all the required derivatives directly from the 
observational data. 

Numerical simulation experiments show that 
ml 1 1 i magni tude accuracy can be achieved — roughly an order 
of magnitude improvement over conventional methods — even 
with bands as broad as 5oG Angstroms (full width at half 
maximum). This is similar to the width of the UBV bands. 
However# much closer spacing (about 200 Angstroms) is 
required than the roughly 1000 Angstrom spacing of the UBV 
bands. Thus# the low accuracy of UBV photometry seems 
primarily to be due to Its violation of the sampling 
theorem# as pointed out earlier [Young# 197A]. 

As the uvby 4-color system is even more under samp I ed# 
one would expect even larger aliasing errors to occur in 
It; and# indeed# errors exceeoing a tenth of a magnitude 
are reported oy Manfroid and Sterxen £1987]# even in careful 
worn where many standards are used and the extinction is 
well determined. 
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Append ix 


The extension of the S trbmgren— Kl ng method to higher-order 

terms requires derivatives of both the smoothed stellar 

spectral Irradiance 1 and the atmospheric extinction A# 

s 

as functions of wavelength. Suppose we sample these 
functions at three unequally spaced wavelengths determined 
by the instrumental filters. We need to know the relative 
spaclngs of the samples (i.e.# the photometric passbands)# 
which are required in the Lagrangian interpolation and 
differentiation formulae. 

Let us suppose that the middle sample (band) is displaced 

a fraction f of the separation of the outer two from their 

midpoint. Thus# f may run from -1/2 at the shortest of the 

three wavelengths co +1/2 at the longest# it would oe 0 if 

the middle band were exactly midway between the others, 

(Obviously# we will try to choose filters that make f small.) 

For the sake of generality# let us use x for the 

independent variable and y for the dependent variable# 

rather than wavelength# spectral irradiance# extinction# 

or any other specific quantity. Our three samples are at 

x # x # and x # and the function values are y # y # 

0 12 0 1 

and y . In terms of f# the middle sample is at 
2 

x * - (x ♦ x ) ♦ f (X -x). (Al) 

1 2 U 2 2 0 

Now# if we want the function value somewhere in the interval 
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from x^ to x^p let us similarly specify the interpolation 

position by a parameter g deflneo like f, so that g runs 

from -1/2 at x to +1/2 at x ; g « f at the middle 
O 2 

sample. In terms of f and g* the interpolating polynomial Is 


y( g) 


l£_Z_l.Lj.£_Z i,Z) 

77 +”1717 "" 


a (f - g) (g + 1/2) 

tf -“1727 y 2 * 


^g_+_l/2 ) (g - 1/2) 

77 ♦ 1 / 2 ) (7 -”I/2) y i 

(A?) 


and its derivative is 


< f + 3/2) 

y'(-l/2) « y 

if + i/2) 0 


if ♦ 1/2) 

♦ ————————— y 

if - 1/2) 2 


at x j 

0 


y'if) 


if - 1/2) 

+ y + 

if + 1/2) 0 


if ♦ 1/2) 


If - 1/2) 


at x^i and 


y‘i+1/2) - - 


if - 1/2) 

77 ”7”I7IT y o * 


if + 1/2) 

if - 1/2) 

2 

f 

<f + 1/2) 

if - 1/2) 

1 


if ♦ 1/2) 

if - 1/2) 


( A3 a ) 


i A3b ) 


if - 3/2) 

♦ y 

(f - 1/2) 2 


( A3c) 


® ^ * • At ail three points* the second derivative Is 
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2 


y" « ♦ --------- -- + T/ - T y f _ £717 y i 


+ — y t 

(f - 1/c) 2 


(A4) 


Independent of g# as the second derivative of a parabola 
Is a constant. 

Now we can write Eq.(ll) In the text as 



2 M(z)A( A ) 
1.066 



a" ♦ ( a * > ♦ 


wher e 


d I n I 
d I n X 


a * * 


d I n A 

d In X 


(A5) 



m" * 

(d In \ ) (o »n \ ) 

and all expressions are evaluated at \ • 

o 

We now set a * In A # for i * 1# 2# and 3 • Then# for 
» i 

g - f (the middle band)/ eq. <a3d) gives 


i» « 


r. (^ z \ 

■r 1 

r tf 

- 1/2) 

2 f 

I- hr) 

J 1 

Lit 

»————— a ♦ 
♦ 1/2) 0 

( f ♦ 1/2) ( f - 1/2) '1 


(f ♦ 1/2) 

77"""I7I7 a 2 


] 


(A6) 
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and there nIIi be similar expressions# based on tqs. (A3a) 
and (A3c)# for the other two bands. For all three# we find 
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Because magnitudes are negative logs to the base 2.512... 

Instead of natural logarithms# the equations for m* and m" 

are similar to these# but contain additional factors of 

~1 . Od 57. . . . The magnitudes# uni Ike the extinction 

coefficients (which are measured on an absolute scale)# 

contain additive zero-point terms due to the instrumental 

sensitivity differences among bands. Thus# in terms of the 

extra-atmospheric monochromatic magnitudes m # m # and 
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for the middle band at g * f # and corresponding equations# 
mutatis mutandis# for the other two bands. Finally# 
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1.0857, 


(A9) 


If- 1/2) 


m 


J \1.0857^ 


at ail three bands. Note tnat the factor of -1.0857... is 

squared here* and hence positive. Also# the magnitudes 

to m here implicitly contain zero-point terms that must 
2 

be evaluated; a constraint such as z ♦ z ♦ z * 0 

0 12 

must be imposed to prevent the matrix of the normal 


equations from oeing singular. 

Thus* the equation of condition for photometric 
reductions is the result of combining Eqs. ( A5 - A9) with 
Eq. (10) in the text. In this combination# note that the 
terms containing a'# a M # m*» and m" always involve either 
second derivatives alone# or squares or products of first 


derivatives# so that the factor 


r ^_/ 1.*** ] 2 

L V" • v v J 


( A10) 


may be removed from all terms. The quantity W in Eq. (A10) 
Plays a role similar to that of the old parameter of the 
same name in cqs. (3.1.56) and (3.1.57) of Young C1974]# but 
the new equations are more exact and Involve fewer 
approximations. In particular# I now evaluate the second 
derivatives of A and I explicitly from the data# as well as 
Keep i ru the lojarlthm function intact# 


The greater complexity of these equations is not a serious 
obstacle to photometric reductions. If we are to approach 
the precision that has long been the prerogative of the 
astr omet r is ts; it is only reasonable that our equations 
must begin to approach theirs in complexity. In any case; 
computers are now so large and fast that there Is no 
difficulty in solving for a slightly larger number of 
parameters! here Me have 3 magnitudes for each star; 3 
extinction coefficients for each night; two independent 
instrumental zero points (which should remain fixed If 
the instrument is well designed and constructed); the three 
bandwidth parameters W; one for each band; and the parameter 
f that specifies the relative band spacing. 

Thus; only the last 6 parameters describe the 
instrumental system. As Manfroid and Heck CI983; 19843 have 
shown; even more Instrumental parameters can be well 
determined if data from several nights are combined. And; 
in fact; numerical experiments show that tnese parameters 
can be determined adequately with a modest number of 
obser vat ions. A particular advantage of this more precise 
model is that the data are represented more closely than 
with the older approximations; so that (If the observations 
are carefully done) the residuals from the least-squares 
fit are smaller. This means that fewer observations per 
parameter are required to reach a given level of precision. 
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However# one must bear in mind that the model is useful 
only if the bands overlap enough to satisfy the sampling 
theor em. 1 have pointed out before (Young# 19743 that no 
existing system does this. Further work is needed to 
determine the necessary spacing; but preliminary numerical 
experiments suggest that bands should be spaced about 1/2 
of their full width at half maximum. Thus# for bands as wide 
as those of the UBV system# a spacing on the order of 200 A is 
sugg e ste d. 
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Table 1. Standard Deviations* Maximum Residuals* and Maximum 
Extrapolation Errors of Extinction Fits. 


Data set 

band 

separ at i on 

Black bodies 

300 

A 

Par abol ic 

300 

A 

Gunn-St ryk er 

100 

A 

Gunn-Stryker 

200 

A 

Gunn-Stryker 

300 

A 

Gunn-Stryxer 

500 

A 


o v er a 1 1 
st d. dev. 

0? 000023 

0.000078 

0.000071 

0.000092 

0. 00020 

0.00067 


m i dd I e band 
max. res i d. 

0*000033 

0.00011 

0.00023 

0.00022 

0.00025 

0.00050 


middle ban 
max. error 

0*0002 

0.0006 

0.0005 

0.0006 

0.0008 

0.0018 
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